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Abstract 

We provide a "first principles" description of scattering from open quantum 
systems subject to a Lindblad-type dynamics. In particular we consider the case 
that the duration of the scattering process is of similar order as the decoherence 
time of the scatterer. Under rather general conditions, the derivations lead to the 
the following new result: The irreversible time-evolution may cause a reduction of 
the system's transition rate being effectuated by scattering. This is tantamount to a 
shortfall of scattering intensity. The possible connection with striking experimen- 
tal results of neutron and electron Compton scattering from protons in condensed 
matter is mentioned. 

Keywords: irreversible dynamics, entanglement, decoherence, neutron Compton scat- 
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1 Introduction 

The counter-intuitive phenomenon of entanglement fTl between two or more quantum 
systems has emerged as the most emblematic feature of quantum mechanics. Experi- 
ments investigating entnaglement, however, are mainly focused on collections of few 
simple (two- or three-level) quantum systems thoroughly isolated from their environ- 
ment (e.g., atoms in high-Q cavities and optical lattices). These experimental condi- 
tions are necessary due to the decoherence of entangled states. In short, decoherence 
refers to the suppression of quantum superpositions caused by the environment. By 
contrast, entanglement in condensed and/or molecular matter at ambient conditions is 
usually assumed to be experimentally inaccessible. However, two new scattering tech- 
niques operating in the sub-femtosecond time scale provided results indicating that 
short-lived entangled states may be measurable in condensed matter even at room tem- 
perature mm. 

In this paper we provide a first-principles treatment of scattering from "small" open 
quantum systems in condense-matter environments, in the "time window" of decoher- 
ence of the scattering system. That is, the focus is in "fast" scattering processes with 
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a duration (usually denoted scattering time, r,c) of the order to the scatterer's decoher- 
ence time, Tdec- This may be considered to represent an "extension" of standard scat- 
tering theory — as applied e.g. to neutron physics ||2][3] or electron scattering JU — in 
which the concepts of entanglement and decoherence play essentially no role. The first 
part of the derivations are analogous to the standard (often denoted) "van Hove formal- 
ism" 111; see also the textbook |3 1. Then a reduced open quantum system, i.e. a micro- 
or mesoscopic system characterized by a set of preferred coordinates, is introduced. 
This corresponds to the "small" physical system that scatters a neutron (electron, etc.) 
with a sufficiently large momentum transfer Its dynamics is described by a simple 
Lindblad-type master equation jT) [8] (which, for the sake of simplicity, contains only 
one Lindblad operator, X), thus including explicitly the effect of decoherence into the 
formalism. 

The striking result of the derivation may be summarized as follows: The irreversible 
time-evolution (owing to the Lindblad operator X) may cause a reduction of the tran- 
sition rate of the system (from its initial to its final state). In "experimental" terms, this 
is tantamount to an effective reduction of the system's cross-section density and thus a 
shortfall of scattering intensity. 



2 Scattering in brief 

We assume an N-body Hamiltonian Htotai = Hq + V with an interaction of the form 

V{r)^Xn{r) (1) 

where n(r) is the particle density operator and A is the rest of the interaction (contact 
potential). For example, in the case of neutron scattering from a system consisting of 
N particles with the same scattering length b one may put 

^w = -^E^('--i^^) (2) 

j 

where V is the volume, JUj is the spatial position of the j-th particle, and 

A^^i., (3) 

m 

m being the neutron mass. (For further details about scattering from "bound" and 
"free" particles, see the textbook [3|.) 

In the interaction picture, the Schrodinger equation is now (putting for simplicity 

h=i) 

idt"^ = \ n{v,t)^ (4) 

with the perturbative solution 

^{t) =^{Q)-i\ f n{Y,t')dt'^{0). (5) 

-'0 
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We write the transition probability W{t) between initial states tpi (with probability 
Pi that the scattering system is in the state ijji) and final states i/jf of the scattering 
system to be given by 



(6) 



It should be noted that ipi and i/jf are eigenstates of the unperturbed N-body Hamilto- 
nian Hq |l3]|2l. This allows us to write the transition probability in the form 



W{t) ^ A' 

where 



/ dt' f dt"Y,{^f I n{r,t')pn{r,t") \ ^f), (7) 
Jo Jo f 



by noting that nt(r, t) — n{r, t), since Rj and r are Hermitian operators. 

In an actual scattering experiment from condensed matter, we do not measure the 
cross-section for a process in which the scattering system goes from a specific initial 
state 4'i to another state ipf, both being unobserved states of the many-body system. 
Therefore, one takes an appropriate average over all these states ID |2l, as done in 
Eq. dSl). 

Furthermore, the initial (ko) and final (ki) momenta of an impinging probe particle 
(neutron) may be assumed to be well defined |[3][2l . Introducing the momentum transfer 
q = ko — ki from the probe particle to the scattering system, the Fourier transform of 
the particle density reads 



n{r, t) = j dqn(q, t) exp{f q • r} 



(9) 



(2^) 

where, in the case of neutron scattering, cf. Eq.©, 

n(q,i) = Eexp{~iq.Rj(t)}. (10) 

Since n(r, t) is Hermitian it holds ri^(q, t) — n(— q, t) and one obtains from Eq. (|6]l 



W{t) = A 



t 
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/ dt' ( dt" V(7/;/|n(q,i')pn(-q,mV'/), dD 
Jo Jo f 



At this stage one traditionally assumes that the sum over V'/ runs over all possible 
eigenstates of Hq which constitute a complete set, i.e. E/|^/;/)('0/| = 1; see Ell]. 
Hence 

Y,{i^fHci,t') pn{-q,t")\^jf) =Tr[n{q,t') pn{~q,t")] (12) 
/ 

where Tr[...] denotes the trace operation. As done in standard theory |l3]|2l, in Eq. ( fT2] l 
one first sums over all final states, keeping the initial state ij^i fixed, and then averages 
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over all tpi (see e.g. f3l, p. 19). The right-hand-side of Eq. (fT2] l contains the density 
operator p of the system before collision, Eq. (O, which is a well known result. 

By introducing a measurement time (the so-called scattering time) Tsc, that is the 
duration of the scattering process, we find 

WiTsc) = /" dt' f dt"Tr[n{q,t')pn{-q,t")] 
Jo JO 

- X'^Tsc f dTTr[n{q,t')pn{-q,t' + t)], (13) 
Jo 

where the stationary property of the correlation function has been used Now one 
can introduce the transition rate, say W, which is defined as 

W = ^^^'"^ = X'' dTTr{n{k,t')pn{--k,t' + t)) 

Tsc Jo 

= dTC{q,T). (14) 

^0 

Here the correlation function 

C(q,t) =rr[n(q,0)pn(-q,i)] (15) 

is introduced, which is analogous to the so-called intermediate function of neutron 
scattering theory Q. 

3 Irreversible dynamics 

We now introduce a set of preferred coordinates { |^)}, cf. [9, JOJ. These are the 
relevant degrees of freedom coupled to the neutron probe. The density matrix needed 
in (flJl i is then the reduced one in the space spanned by these states, and it is obtained 
by tracing out the (huge number of the) remaining degrees of freedom belonging to 
the "environment" of the microscopic scattering system (e.g. a proton and its adjacent 
particles). To simplify notations, we denote this reduced density matrix by p too. 

In the subspace spanned by the preferred coordinates (also denoted 'pointer basis'), 
we assume the relevant density matrix to obey a Lindblad-type equation of the form 

mm 

dtp=-i[H,p]+np = Cp (16) 

with the formal solution 

Pit) = e^V(O). (17) 
Let us look at a time-dependent expectation value 

{Ait)) = Tr ipit)A) = Tr (e^V(O)A) = Tr (p(0)e^'*A) , (18) 

where we define by setting 

TriiCX)Y)=TriXiC^Y)) . (19) 
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Thus we obtain a Lindblad time evolution for the operators too by writing 

dtA{t)^C^A{t). (20) 

This form was actually the original Lindblad result. Note that this works as long as C 
does not depend on time. For time-dependent generators of the evolution, a somewhat 
more elaborate scheme is needed. 

Now we find that we may use this formalism to calculate correlation functions like 
the one in (fT4l i. We write 

{A{t)B) = Tr [p{Q) (e^''A) b\ = Tr [Ae^* (i?p(0))] = Tr {ApB{t)) , (21) 

where pB{t), as defined in Eqs. (l2Tl i. obeys the equation 

dtPBit) = Cpsit) (22) 

and the initial condition 

PsiO) = Bp{Q). (23) 

Thus, except for the initial condition, we have to solve the same equation of motion as 
for the density matrix ( fTSI l. 

4 Application to scattering 

We here assume a simple Lindblad-type ansatz for the master equation in the relevant 
subspace. We set 

dtp = -I [H, p] - K [X, [X, p]] - Cp, (24) 

where K > Q and H is the reduced (or relevant Hamiltonian) of a microscopic or 
mesoscopic scattering system and the double commutator term describes decoherence. 
For simplicity of the further calculations, we here assume that 



(25) 



With Eq. (flST l we have 

C(q, t) ^ Tr[n(q, 0) p ri(-q, r)] = Tr (^n(q, 0) pe^^'*n(-q, 0) 

- Tr(n(-q,0)e^*(n(q,0)p)) . (26) 
This is equivalent with the expression 

C(q,r) = I n(-q,0) | I Pn{t) \ 0- (27) 



With the equation (l241 i. one easily finds the well known solution 

{^'\Pn{t)\0 = expH(% -^«)i]exp[-A-(e'-0'i 

X \n{q,0)pi0)\O- (28) 
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Inserting this into the expression (fT4] l for the transition rate we find 

x(e|n(-q,0)U')(e'hKq,0)p(0)U)- (29) 

5 Decoherence and decrease of transition rate 

Obviously, the decoherence-free limit of this result, i.e. with K = 0, corresponds to 
the conventional result of scattering theory. 

The oscillating factors exp [—i {£^> — £^)t] are characteristic for the 'unitary-type' 
dynamics caused by the commutator part —i [H, p] of the master equation ( |24] | for the 
reduced (or: relevant) density matrix p. These factors have the absolute value 1 and do 
not affect the numerical value of the transition rate. 

On the other hand, the restrictive factors exp(— if — t) < 1, which are due 
to the decoherence, can be seen to cause a decrease of the transition rate and thus of 
the associated cross-section. This can be illustrated in physical terms as follows: 

Let us first assume that the reduced density operator p(0) can be chosen to be 
diagonal in the preferred (^—representation (which corresponds to the usual random 
phase approximation at t = 0). Then each term of Eq. ( |29l l is of the form 

(^|n(-q,0)|e')(e>(q,0)p(0)|O = 

(^K-q,o)U')(^'Kq,o)IO(C|p(o)IO = 

|(e|n(-q,0)|OP(^|p(0)IO>0. (30) 

The last inequality is valid because it holds (C|p(0)|O > 0. 

If the assumed diagonal form of p(0) would be considered as being 'too strong', 
one may note the following: The exponentials exp(— A' {£,' ~ ^) t) due to decoherence 
imply that only terms with ^ ~ ^' contribute significantly to the transition rate. Thus 
we may conclude that, by continuity, all associated terms with ^ sa ^' in Eq. ( |29] | should 
be positive, too. The further terms with ^ being much different from ^' can be positive 
or negative. But they may be approximately neglected, since they decay very fast and 
thus contribute less significantly to W. 

The main conclusion from the preceding considerations is that the time average 
in Eq. (|29] | always decreases the numerical value of W = W{tsc)/tsc, due to the 
presence of the exponential factors exp{~K — ^) t) < 1. In other words, the 
effect of decoherence during the experimental time window t^c plays a crucial role 
in the scattering process and leads to an 'anomalous' decrease of the transition rate 
and the associated scattering intensity. This result is in line with that of Ref. ifTTl . 
which investigated the standard expression of the double differential cross-section of 
neutron scattering theory by ad hoc assuming decoherence of final and initial states of 
the scatterer. 

Very interesting is also the conclusion that, in the limit of very slow decoherence 
(K 0), this 'anomaly' disappears, i.e. the scattering results are expected to agree 



W = (ir ^exp[-i(% -f^)^] exp 
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with conventional theoretical expectations. This is contrary to the associated prediction 
of the theoretical model of Ref. Iil2il 

6 Additional remarks 

A related effect (i.e., a shortfall of scattering intensity) was observed in recent neutron- 
proton Compton scattering (NCS) and electron-proton Compton scattering (ECS) ex- 
periments in condensed matter ||5]|6l, in which the experimental scattering time lies in 
the sub-femtosecond time scale. This coincides with the characteristic time of elec- 
tronic re-arrangements accompanying the breaking (or formation) of a chemical bond. 
Note that in these experiments the energy transferred to a proton is large enough to 
break the bond (C-H and O-H). 

Some remarks about the possible selection, definition and/or physical meaning of 
the preferred coordinates may be appropriate. In the case of conventional NCS theory, 
for example, one uses momentum eigenstates of the scattering particle (e.g. proton) — 
as well as for the neutron — as the appropriate basis 1 13]. In the light of the preceding 
derivations, however, one may observe the following: Due to the strong (Coulomb) 
interactions of the scattering proton with its adjacent particles (electrons, and probably 
also other nuclei), { | ^ } can not be one-body states but they should rather be considered 
to represent momentum states being strongly "dressed" (and entangled) with degrees 
of freedom of adjacent particles. 

Further work will deal with the more general — and experimentally relevant — 
case, in which the preferred states {|^)} are not eigenstates of the "reduced" energy 
and Lindbald operators, Eqs. (|25l l. In that case, the result of Eq. (|29] l will become less 
simple. 
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